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We report on the calculation of the next-to-leading order QCD corrections to the production of 
W-boson pairs in association with a hard jet at the Tevatron and the LHC, which is an important 
source of background for Higgs and new-physics searches. Leptonic decays of the W bosons are 
included by applying an improved version of the narrow-width approximation that treats the W 
bosons as on-shell particles, but keeps the information on the W spin. A selection of differential 
NLO QCD cross sections is provided both for the LHC and the Tevatron. The QCD corrections 
stabilize the LO prediction for the cross section with respect to scale variations. The differential 
LO cross sections are generally not simply rescaled by the corrections. Their shapes are particu- 
larly distorted if an additional energy scale is involved. 
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1. Introduction 

The search for new-physics particles — including the Standard Model Higgs boson — will be 
the primary task in high-energy physics in the era of the LHC. The extremely complicated hadron 
collider environment does not only require sufficiently precise predictions for new-physics signals, 
but also for many complicated background reactions that cannot entirely be measured from data. 
Among such background processes, several involve three, four, or even more particles in the final 
state, rendering the necessary next-to-leading-order (NLO) calculations in QCD very complicated. 
This problem lead to the creation of an "experimenters' wishlist for NLO calculations" [[I], ^ ||] 
that were still missing at that time, but are required for successful LHC analyses. The process 
pp — )• W+W^+jet + X made it to the top of this list. Meanwhile the 2^3 particle processes and 
also some of the 2 — > 4 particle processes [Q, ^ on the list have been evaluated at NLO QCD. 
Moreover, benchmark results for the virtual corrections have been presented for a specific phase- 
space point for all 2 — > 4 processes on the list in Ref. 

The process of WW-f-jet production is an important source for background to the production 
of a Higgs boson that subsequently decays into a W-boson pair, where additional jet activity might 
arise from the production. WW-i-jet production delivers also potential background to new-physics 
searches, such as supersymmetric particles, because of leptons and missing transverse momentum 
from the W decays. Besides the process is interesting in its own right, since W-pair production 
processes enable a direct analysis of the non-abelian gauge-boson self-interactions, and a large 
fraction of W pairs will show up with additional jet activity at the LHC. Last but not least WW-f-jet 
at NLO also delivers the real-virtual contributions to the next-to-next-to-leading-order (NNLO) 
calculation of W-pair production, for which further building blocks are presented in Ref. [Q]. 

Here we report on the calculation of the process pp/pp W+W^+jet + X in NLO QCD in- 
cluding leptonic W-boson decays. Results of this calculation have been published in Refs. 



Parallel to our work, another NLO study [10] of pp ^ W+W +jet + X at the LHC appeared 



Moreover, NLO QCD corrections to the related processes pp^Wy+jet + X [ p| ] and 



pp — ZZ+jet + X 1 12] have been calculated recently. 



2. Details of the NLO calculation 

At leading order (LO), hadronic WW+jet production receives contributions from the partonic 
processes qq — )• W+W^g, qg — > W+W^q, and gq — )• W+W^q, where q stands for up- or down- 
type quarks. Note that the amplitudes for q = u, d are not the same, even for vanishing light-quark 
masses. All three channels are related by crossing symmetry. 

The leptonic W decays are implemented by means of an improved narrow-width approxima- 
tion (NWA) that treats the W bosons as on-shell particles, but keeps the spin correlations between 
production and decay processes. In this way, a significantly better approximation of the full calcu- 
lation is achieved, which can be read off the comparison of the sample LO distributions in Figure 

In order to prove the correctness of our results we have evaluated each ingredient twice us- 
ing independent calculations based — as far as possible — on different methods, yielding results in 
mutual agreement. 
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Figure 1: Comparison of W-decay descriptions in the distributions of the transverse momentum of e+ (left 
plot), the pseudo-rapidity of e+ (central plot), and the azimuthal angle between the two decay leptons (right 
plot). The LO cross sections are evaluated at = /Xfact = Mien = for the full calculation, the naive NWA, 
and the improved NWA. (Taken from Ref. [^.) 



2.1 Virtual corrections 



Version 1 of the virtual corrections is essentially obtained as for the related processes of ttH 



[13] and tt+jet [14| production. The Feynman diagrams are generated with FeynArts 1.0 [15] and 
further processed with in-house Mathematica routines, which automatically create an output in For- 
tran. The IR divergences (soft and collinear) are analytically separated from the finite remainder in 
terms of triangle subdiagrams, as described in Refs. [Oj, 1^. This separation, in particular, allows 
for a transparent evaluation of so-called rational terms that originate from D-dependent terms multi- 
plying IR divergences, which appear as single or double poles in £. As generally shown in Ref. 
after properly separating IR from UV divergences such rational terms originating from IR diver- 
gences completely cancel; this general result is confirmed in our explicit calculation. For the results 
presented in Ref. the pentagon tensor integrals were directly reduced to box integrals following 

to scalar 



Ref. [jlTp , while box and lower-point integrals were reduced a la Passarino-Veltman 
integrals. This procedure completely avoids inverse Gram determinants of external momenta in 
the reduction step from 5-point to 4-point integrals, but the reduction of box and lower-point ten- 
sor integrals involves such inverse determinants via the Passarino-Veltman algorithm. Although 
these inverse determinants jeopardize the numerical stability in regions where such determinants 
are small, sufficient numerical stability was already achieved. Meanwhile the tensor reduction has 



been further improved using the methods of Ref. []19|]. The scalar one-loop integrals are either 
calculated analytically or using the results of Refs. [||, ||, H). 

Version 2 of the evaluation of loop diagrams starts with the generation of diagrams and am- 
plitudes via FeynArts 3.4 [23] which are then further manipulated with FormCalc 6.0 [|^] and 
eventually automatically translated into Fortran code. The whole reduction of tensor to scalar inte- 
grals is done with the help of the LoopTools library [24|, which employs the method of Ref. (IF 



for the 5-point tensor integrals, Passarino-Veltman [18] reduction for the lower-point tensors, and 
the EFpackage [^, ^] for the evaluation of regular scalar integrals. The dimensionally regularized 
soft or collinear singular 3- and 4-point integrals had to be added to this library. To this end, the 
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Figure 2: Scale dependence of the WW+jet cross sections with W decays included and further cuts applied 
according to Ref. [^. For the LHC setup, the results are given for ;?T,jet > lOOGeV (left plot). For the 
Tevatron we show results for PTjet > 50GeV (right plot). 



explicit results of Ref. [16] for the vertex and of Ref. [27] for the box integrals (with appropriate 
analytical continuations) are taken. 

2.2 Real corrections 

The matrix elements for the real corrections are given by the processes W+W^qqgg and 
W+W^qqq'q' with a large variety of flavour insertions for the light quarks q and q'. The 
partonic processes are obtained from these matrix elements by all possible crossings of quarks and 
gluons into the initial state. The evaluation of the real-emission amplitudes is performed in two 
independent ways. In one approach we apply the Weyl-van-der-Waerden formalism (as described 
in Ref. [E^]). The other one is based on Madgraph [E^ generated code. Both evaluations employ 



(independent implementations of) the dipole subtraction formalism [ pO| ] for the extraction of IR 
singularities and for their combination with the virtual corrections. 

In one calculation the phase-space integration is performed by a multi-channel Monte Carlo 



integrator [31] with weight optimization [32] written in C++, which is constructed similar to 



RacoonWW [33, 34]. The second calculation uses a simple mapping where the phase space is 
generated from a sequential splitting. 



3. Numerical results 



We consistently use the CTEQ6 [ pSQ set of parton distribution functions (PDFs), i.e. we take 
CTEQ6L1 PDFs with a 1-loop running in LO and CTEQ6M PDFs with a 2-loop running in 
NLO. The complete setup we used for our numerical calculations is precisely defined in Ref. [^], 
where a large variety of additional results is provided. 

Figure ^ shows the scale dependence of the NLO cross section for the LHC and the Tevatron. 
The QCD corrections stabilize the LO prediction for the WW-f-jet cross section considerably with 
respect to a variation of the factorization and renormalization scales which we identify with each 
other. At the LHC, this stabilization of the prediction, however, requires a veto on a second hard jet. 
Otherwise the production of final states with WW-H2jets, which yields a LO component of the NLO 
correction, introduces again a large scale dependence. In Figure ^ a sample of NLO distributions is 
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Figure 3: Differential cross sections for WW+jet with decays included in the improved NWA at the LHC: 
The LO and NLO distributions are shown for ji = jUfact = Mren =Myj. The distributions of the transverse 
momentum of e+ (left plot), the pseudo-rapidity of e+ (central plot), and the azimuthal angle between the 
two decay leptons (right plot) are depicted. The bands in the /T-factors refer to a variation of by a factor 
of 2 in the NLO quantities. (Taken from Ref. [^.) 



provided for the LHC setup. At the LHC the psuedo-rapidity distributions in the dominant region 
and also the distributions in the angles between the two charged leptons have an almost constant 
^-factor of about 1.3 (inclusive cross-section definition); for the exclusive cross-section definition 
the corrections are even smaller and rather close to L The pj spectra, on the other hand, show 
a much more phase-space-dependent /T-factor with the exclusive cross-section definition showing 
an even larger dependence than the inclusive one. This is not surprising since the pj introduces 
an additional scale which could introduce potentially large logarithms which are badly treated by 
a constant renormalization scale. At the Tevatron our findings are similar [^, ^]. We note that 
the almost constant /T-factor which holds for a remarkable number of distributions has also been 



observed in Ref. [10]. 
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